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Abstract This contribution is not about the quality of the 
agreement between stellar models computed by CESAM and 
CLES codes, but more interesting, on what ESTA-Task 1 run 
has taught us about these codes and about the input physics 
they use. We also quantify the effects of different implemen- 
tations of the same physics on the seismic properties of the 
stellar models, that in fact is the main aim of ESTA experi- 
ments. 
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1 Introduction 

The goal of ESTA-Task 1 experiment is to check the evolu- 
tion codes and, if necessary, to impro ve them. The results of 
Tas k 1 comparisons were pr esented in lMonteiro et al\ (2006) 
and iLebreton et al.\ d2008al) for a set of stellar models rep- 
resentative of potential CoRoT targets. The models calcu- 
lated for TASK 1 were based on rather simple input physics. 
Moreover, a great effort was done to reduce at maximum 
the differences between computations by fixing the values 
of fundamenta l consta nts and the physics to be used (see 
ILebreton et all l2008bl) . In spite of that, some differences 
among stellar models computed by different codes persist. 

In CESAM and CLES computations we paid attention to 
adopt, not only the same fundamental constants and metal 
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mixture dGrevesse and Noefsl I1993L thereafter GN93), but 
also the same isotopic ratios and atomic mass values. Nev- 
ertheless, even if the same metal mixture, opacity tables and 
equation of state were adopted, there is still some freedom 
on their implementation. In this paper we analyze these dif- 
ferent implementations and estimate the consequent effects 
on the stellar structure and on the seismic properties of the 
theoretical models. In section |2] we study the equation of 
state and in section [3] the differences in the opacity tables. 
The nuclear reaction rates are discussed in Sect. [4] and the 
effect of different surface boundary conditions in Sect.|5] Fi- 
nally, in Sect. [6] we analyze the differences due to different 
numerical techniques. 



2 Equation of State 



As fix ed in ESTA we used the OPAL2001 dRogers and NavfonovL 
2002) equation of state which is provided in a tabular form. 
In CESAM the quantities: density, p, internal energy, E, the 
compressibilities %T = {dlnP/dlnT) p and^p = (dlnP/ dlnp)j, 
the adiabatic indices il, li/C^ - 1), I3 — 1, and the specific 
heat at constant volume Cy, are obtained from the variables 
P, T, X and Z (respectively pressure, temperature, hydro- 
gen and heavy element mass fraction) using the interpola- 
tion package provided on the OPAL web site, and the spe- 
cific heat at constant pressure (C p ) is derived from (I3 — 1). 
On the other hand, CLES interpolates only Cy, P, % p and %T 
in the OPAL EOS tables by a method ensuring the continuity 
of first derivatives at cell boundaries in the four-dimensional 
space defined by the variables p,T,X and Z. The other ther- 
modynamic quantities Fi , (F3 — 1) and C p are derived from 
the values of Cy, P, % p and %T by means of the thermody- 
namic relations. 

As a first step we want to disentangle the differences in 
the thermodynamic quantities from their effects on the stel- 
lar structure. We estimate therefore the intrinsic differences 
between the equation of state used in CESAM and in CLES. 
To this purpose we computed the differences between the 
thermodynamic quantities from the corresponding EoS rou- 
tines, for a stellar structure defined by p, T, X and Z values. 
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Fig. 1 Relative differences in the F values provided for a given (p,T) 
structure by CESAM and CLES EOS routines. The thick lines corre- 
spond to Fi values derived from Cy tabulated values while the thin ones 
are derived from the tabulated F . Solid lines corresponds to a 2 Mr... 
model with a central hydrogen mass fraction X c = 0.50, and dashed 
ones to a 0.9 M@ star in the middle of the main-sequence. 
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Fig. 2 Frequency differences between CESAM and CLES 2 M mod- 
els. The two different curves correspond to CLES models computed by 
using the two different EoS tables (see text). 
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In Fig. Q] (thick lines) we show the result of Fi comparison 
for two different stellar models, a 2 M model with a mass 
fraction of hydrogen in the center X c = 0.50 (solid line) and 
a 0.9 Mq model with X c = 0.35 (dashed line). By comparing 
also the other thermodynamic quantities we found that the 
largest discrepancies between CESAM and CLES EoS occurs 
for logF < 5 (corresponding to the partial He and H ion- 
ization regions), and they are, at maximum, of the order of 
2% for r u and of 5% for V ad and C p . By using the OPAL 
interpolation routine in CLES, we verified that the different 
interpolation schemes used in CESAM and CLES can only ac- 
count for an uncertainty of 0.05% in P, 0.2% inFi, and 0.5% 
in an d C p . These remaining differences are probably ex- 
plained by the fact that CESAM uses as variables (P, T) and 
uses subroutine rhoofp of OPAL-package to transform (P, T) 
into (p,F), while CLES uses directly (p,F). Nevertheless, 
those discrepancies are an order of magnitude smaller than 
the differences between CESAM and CLES EoS. 

As it was already noted by Boothro vd and Sackmannl (120031) ., . . 
some inconsistencies existed between thermodyna mic quan- th e internal emr ^ as tabulat ed in OPAL EOS, showed that 
tities tabulated in OPAL EOS: "for the OPAL EOS \R02ers et alf & 0PAL tabulated Cy was affected by a large inaccuracy. 
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Fig. 3 As Fig.|2]for 0.9 M e model. 



1996), we found that there were significant inconsistencies 
when we compared their tabulated values ofFy, Fx/ {pi — 1), 
and (I3 — 1) to the values calculated from their tabulated 
values of P, Cy, Xp> an d Xt- ■■■ Preliminary tests indicate 
that this OPAL2001 EOS has larger but smoother inconsis- 
tencies in its tabulated thermodynamic quantities...". As a 
consequence of these inconsistencies, the choice of the ba- 
sic thermodynamic quantities is not irrelevant, and it was 
shown by Roxburgh (2005, priv. communication), that the 
choice done in CLES was the worst one. A direct compari- 
son with the values of Cy computed from the derivative of 



The OPAL team acknowledged afterwards the Cy -issue 
and recommended not to use it. The EoS tables used in CLES 
have then been changed by replacing the tabulated Cy value 
by that obtained from the tabulated values of P, % p , %j, 
and fj. The remaining discrepancies (~ 0.2%) between the 
CESAM Fi values and those from the new CLES-EoS table 
(hereafter called CLES-EoS-F to tell it apart from the origi- 
nal one CLES-EoS-Cy) are due to the different interpolation 
routine. As shown in Fig.[T](thin lines) the discrepancies are 
much smaller for a solar like than for a 2 Mq model and they 
appear mainly in the ionization regions. 
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Concerning the quantities that in CLES are obtained from 
thermodynamic relations and in CESAM from interpolation 
in OPAL tables, the differences come in part from the inter- 
polation routine and in part from the remaining, even if much 
smaller, inconsistencies between the tabulated values of Fi , 
rz/{r% — 1) = V^ 1 , and (Ij — 1). For instance, the values of 
Cy derived from ti may differ by 0.5% from the correspond- 
ing value obtained from (F3 — 1), and that occurs always in 
the H and He ionization regions. The problem is that even if 
the thermodynamic relations to derive the adiabatic indices 
seem more physical, there is some numerical incoherence. In 
fact, the derivatives of interpolated (very often polynomial) 
quantities do not fit in those of the interpolated functions 
(whose behavior is far from polynomial one). 

All t he CLES models involved in Task 1 and Task 3 com- 
parisons (Lebr eton et a/.Ll2008al : lMonteiro et fl/.Ll2008h were 
recomputed with CLES -E0S-F1, but the model s used for com- 
parisons presented in iMonteiro et al\ d2006h were not. In 
fact, most of the frequency differences found in that paper 
came from CLES-EoS-Cv . The effect of EoS differences on 
the seismic properties are illustrated in Fig. [2] for a 2 A/© 
model and in Fig. [3] for the solar like model. In those fig- 
ures we plot the frequency differences of I = modes for 
CESAM models and two types of CLES ones: those computed 
with EoS-Cv (dashed lines) and those computed with EoS-/] 
(solid line). The period of the oscillatory signature shown 
by Av (Vcles — Vcesam) m Fig. [2] is related to the acous- 
tic depth where models differ. A Fourier transform of Av 
shows clearly that the oscillation is linked to the F differ- 
en ces. Moreover the co mparison of Fig. [3] with the Fig.6 
in lMonteiro et al\ (120061) confirms that also for the Case 1.1 
model (see e.g. iLebreton et al.l [2008a), the maximum dif- 
ference of almost 2 jUHz betw een CESAM and CLES models 
found by Mo nteiro et at l (l2006h was due to the inconsistency 
between the tabulated Cy and adiabatic indices. 



3 Opacities 

ESTA specifications require the use of OPAL96 opacity ta- 
bles dlglesias and Rogerslll996l) complemented at low tem- 
peratures by the Alexander and Ferguson ( 1994) (thereafter 
AF94) tables. CESAM uses OPAL tables provided by C. Igle- 
sias, prior to their availability on the web site, and inter- 
polates in the opacity tables by means of a four-point La- 
grangian interpolation. The OPAL opacity tables used by 
CLES were picked up later on the OPAL web site and smoothed 
according to the prescription found there (xztrin21.f routine), 
we will call them thereafter OPAL96-S. Furthermore, the 
interpolation method in CLES opacity routine is the same as 
that used in EoS table interpolation. In both codes the metal 
mixture adopted in the opacity tables is the GN93 one. 

To disentangle the differences in the opacity computa- 
tions from the differences in the stellar structure, we pro- 
ceed as in EoS table analysis, that is, we estimate the intrin- 
sic differences in the opacity (jc) by comparing the K values 
provided by CESAM and by CLES routines for the same stel- 
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Fig. 4 Relative differences in the opacity values provided for a given 
(p,T) structure by CESAM and CLES opacity routines. The thick lines 
correspond to the jc values derived from the smoothed OPAL tables, 
while the thin ones were obtained by including in the CLES opacity 
routine the OPAL tables without smoothing. Solid lines corresponds to 
a 2 M Q model with X c = 0.50, and dashed ones to a 0.9 M Q star in the 
middle of the main-sequence. 
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Fig. 5 The same differences as in Fig. [4] but plotted as a function of 
the relative stellar radius. 



lar structure. The results of these comparisons are shown in 
Figs. [4] and |5j where we plot for two different stellar mod- 
els the opacity relative differences ( JCcesam ~ Kcles ) / K" as 
a function of the local temperature and of the relative ra- 
dius. From comparisons of different models it results that 
the opacity discrepancies depend on the mass of the stel- 
lar model and, for a given mass, on the evolutionary state 
as well. Moreover, a peaked feature at log T ~ 4 which can 
reach values of the order of 5%, appears in all the com- 
parisons. This is a consequence of the differences between 
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OPAL and AF94 opacities in the domain [9000 K-12000 K] 
and of the different method used in CLES and CESAM to as- 
semble AF9 4 and OPAL tables. C LES uses the procedure 
described in Scuflai re et al. I (120081) that ensures a smooth 
passage between both tables, while CESAM searches for the 
point of minimum discrepancy between OPAL and AF94. 
In the interior regions the differences between CESAM and 
CLES opacities do not present the oscillatory behavior that 
we would expect if these differences resulted from the inter- 
polation schemes. On the contrary, the CESAM opacities are 
systematically larger (by 1-2%) than the CLES ones in the 
region logT G [5.5,7] of 2 M G model. 

Even if the metal mixture to be used in opacity compu- 
tations is fixed (GN93), there may be some uncertainties in 
its definition. For instance, OPAL uses atomic masses that 
do not correspond to the value s given by the isotopic ratios 
in lAnders and Grevessd dl989h . In particular there is a dif- 
ference of 0.5% for Neon, and 10% for Argon. Moreover, 
OPAL opacity tables are computed for 19 elements, while 
the GN93 mixture contains 23 elements. There are two op- 
tions: either to ignore the mass fraction of F, Sc, V, and 
Co, or to allot the abundances of these elements among the 
close neighbors. We have analyzed the effects of these un- 
certainties on the opacity values, but they turned out to be 
of the same order of the accuracy in OPAL data (0.1-0.2%). 
Hence, they cannot account for the discrepancy between CE- 
SAM and CLES opacities. 

The other important difference between CESAM and CLES 
opacities is on whether they use the OPAL smoothing rou- 
tine or not. In fact, the OPAL opacity tables are affected 
by somewhat random numerical errors of a few percent. To 
overcome undesirable effects the OPAL web site suggests to 



pass the original tabular data through a smoothing filter be- 
fore interpolating for Z, X, log T, and R (with R = p /T^). A 
direct comparison between the original and smoothed opac- 
ity values have shown a difference larger than 2% for log/? = 
—4 and —3.5 and logT G [5.5,7]. These differences decrease 
for larger and smaller values of log/?. 

We have computed new opacity tables for the CLES opac- 
ity routine without passing through the smoothing filter (CLES- 
OPAL96 instead of CLES-OPAL96-S). The comparison between 
CESAM and new CLES opacity computations are also shown 
in Figs. |U|5] (thin lines). We note that when both codes use 
similar OPAL96 tables, the discrepancies in the internal re- 
gions almost disappear. The remaining differences are due 
to the interpolation schemes and to the small differences in 
GN93 definition. The feature at log T ~ 4 is still there since 
the method used in CLES to assemble AF94 and OPAL96 
tables is the same as in CLES-OPAL96-S. 

At variance with the EOS tables, where an error was 
detected and acknowledged by the OPAL team, we do not 
have any argument to prefer the smoothed to the original 
opacity tables, and we think that the differences between 
both groups of results must be considered an estimate of the 
precision of current stellar modeling. Therefore, the Liege 
group decided to provide fo r Task 1 and Task 3 compar- 
isons (Lebreton et al ., 2008a) the modes computed with the 
standard tables in CLES, that is OPAL96-S. A part of the 
differenc es between CESAM and CLES models that were re- 
ported in lLebreton et al.\ (I2008al) should be hence due to the 
opacity tables we used. In order to estimate these effects we 
have re-computed with CLES and OPAL96 tables (without 
smoot hing) the models for a ll the cases in TASK 1 (see Ta- 
ble 1. Lebreton et al.l l2008al) . In the next three sections we 
present the effect of the opacity uncertainty on: the global 
parameters, the stellar structure, and on the seismic proper- 
ties. 
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Fig. 6 Evolutionary tracks for stellar parameters corresponding to the 
case CI. 5 in Task 1, and also without overshooting. Solid thick lines: 
CLES models with the default opacity tables (OPAL96-S); solid thin 
lines: CESAM models ; dotted lines: CLES models where the opacity 
tables have been recomputed without using the opal smoothing filter 
(OPAL96). 



3.1 Effects on global stellar parameters 

In general, the change of opacity tables decreases the dis- 
crepancies between the stellar global parameters provided 
by CLES and CESAM. In table 1 we collect the differences 
(in percent) in radius, luminosity, central density, and central 
temperature between pairs of models for the Task 1 targets. 
Columns labeled A give the differences Xcles-opal96-s ~ 
^CLES-OPAL96> ar, d columns B and C the differences with 
CESAM, that is, X C les-opal96-s — ^cesam and X C les-opal96 — 
^CESAM respectively. 

We note that for the most evolved models (CI. 3 and 
CI. 5), the effect on the radius of changing the CLES opacity 
tables is small, and that the agreement with CESAM gets even 
worse than with the original tables. There is however a sig- 
nificant decrease of the luminosity discrepancy. For the cases 
CI .4 and CI .6 (PMS and ZAMS respectively) the change from 
OPAL96-S to OPAL96 is particularly effective, leading to 
a decrease of AR and AL by a factor 4 and 6 respectively 
for CI. 4, and dropping the discrepancy to values lower than 
0.003% for CI. 6. 
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Fig. 7 Effect of opacity uncertainties on the stellar radius and luminos- 
ity along the main-sequence evolution of a 2 M@ model. Thin lines cor- 
respond to differences between CESAM and CLES-OPAL96-S global 
parameters, and thick ones to differences between CESAM and CLES- 
OPAL96 ones. The dotted lines refer to differences in radius (lower 
curve) and in luminosity (upper curve) between CLES models com- 
puted with OPAL96-S and OPAL96 opacity tables. 



We also studied the effect of opacity tables on the main- 
sequence evolution of a 2 M© star (parameters correspond- 
ing to C1.5). As shown in Fig.|6]ihe HR location of the CLES- 
OPAL96-S evolutionary track is significantly modified by 
adopting OPAL96 tables, and except for the second gravi- 
tational contraction, the new track coincides quite well with 
the CESAM one. The discrepancies in radius and luminosity 
along the MS as well as the effect of opacity tables on their 
values are shown in Fig. [7] These discrepancies are signif- 
icantly reduced by switching from OPAL96-S to OPAL96, 
nevertheless the difference in the stellar radius at the end of 
MS phase (X c < 0.2) is unchanged. That is not due to the 
uncertainties in opacity as shown by the dotted lines that 
correspond to the differences between CLES models com- 
puted with the two different opacity tables. Neither it is a 
consequence of the treatment of overshooting since stellar 
models computed without overshooting for the same stel- 
lar parameters show similar discrepancies. The reason is in 
the treatment of the borders of convective regions in CLES 
that l eads to a sort of "numerical diffusion" dScuflaire et all 
120081) . For the stellar evolution, this diffusion at the border of 
the convective core works as a slightly larger overshooting. 
A significant increase of the number of mesh points used for 
computing the models reduces the numerical diffusion and 
improves the agreement CLES-CESAM. 



3.2 Effects on the stellar structure 



To anal yze to which extent the differences reported in lLebreton 
(2008a) come from the uncertainty in the opacity, we com- 



puted for each TASK 1 model the local differences in the 
physical variables at fixed relative mass and at fixed rela- 
tive radius. To this purpose we use the so-called diff-fgong.d 
routine in the ADIPLS packageQ. In Figs. [8] and [5] we plot 
the logarithmic differences of the sound speed, c, and pres- 
sure, P, for the stellar interior (two left panels) and of sound 
speed and the adiabatic exponent F in the external layers. In 
each panel there are three different curves that correspond 
to the comparisons labeled A, B, and C in previous s ection. 
So, if the differences shown in lLebreton et all (2008a) come 
from the differences in the opacity tables, the solid and dot- 
ted lines should be close to each other. We call again the at- 
tention to the improvement got for the cases C1.4 and CI. 6. 



3.3 Effects on the frequencies 

Since the frequency of p-modes firstly depends on the stellar 
radius, the improved agreement between CESAM and CLES 
models that we obtained by changing from OPAL96-S to 
OPAL96 implies also a decrease in the frequency discrepan- 
cies. The values of Av (Vcles — Vcesam) change from 7.25 
to 6.35 /iHz (at 5300 /iHz) for the case Cl.l; from 7 to 
4.5 /iHz (at 4000 /iHz) for CI. 2; from -5 to -4.5 /iHz at 
1000 /iHz (and from -11 to -10 at 2300 /iHz) for C1.3; from 
3.7 to -1 /iHz at 2200 /iHz for C1.4. The improvement is 
only of 0.3 /iHz for the case C1.5, and Av at 1200 /iHz is 
of the order of 7.5 /iHz. For CI. 6 the initial Av ~ 3 /iHz 
at 1500 /iHz drops to values lower than 0.05 /iHz, and for 
CI. 7, Av change from 1.6 /iHz to 1.2 /iHz (at 600 /iHz). 

By comparing frequencies that have been scaled to the 
same radius we remove the effect of AR and make appear 
the differences due to discrepancies in the stellar structure. 
This was done in Fig.[K)]where we plot the frequency differ- 
ences for £ = and 1 and for the cases considered in Task 1 . 
There are two curves in each panel, one corresponding to 
the difference Vctes-opat .g<s-s - Vce sam (that is, that ap- 
pearing also in Leb reton et all 1200 8a), and the second one 
corresponding to v C les-opal96 - Vcesam- 

The role of opacity on the oscillation frequencies is a in- 
tricate problem since the variations of K lead to changes of 
the temperature structure in the star, and therefore also of 
the value of F- As can be seen in Fig. [9] for the case 1.5, 
the differences in the outer layers might even increase for 
the model computed with similar opacity tables (OPAL96), 
and as consequence, the frequency differences scaled to the 
same radius (Fig. \W[ lower-left panel) show a discrepancy 
even larger than with OPAL96-S. The oscillatory behavior 
is produced by the peak in 51nc at r/R ~ 0.997. A compar- 
ison between CLES-OPAL96 and CLES-OPAL96-S clearly 
shows the same oscillatory behavior, but the absolute fre- 
quency difference is only ~ 0.3/iHz (both models have sim- 
ilar radius). In the same way, CESAM and CLES 2 M Q models 
at X c = 0.50 (whose radius differ by less than 3 x 10~ 4 /?* 
once CLES adopts OPAL96 tables) show frequency differ- 
ences of the order of 0.4 /iHz at 1200 /iHz. After normaliz- 
et all 

1 http://www.corot.pt/ntools 
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Table 1 Global parameter differences, in percent, between Task 1 -target models computed by CLES with OPAL opacity tables in the two ver- 
sions: smoothed (OPAL96-S) and not smoothed (OPAL96). Stellar radius (AR), luminosity (AL), central density (Ap c ) and central temperature 
(AT C ). Columns labeled A give the differences X C les-opal96-s — ^cles-opal96, and columns B and C the differences with CESAM, that is, 
^cles-opal96-s -^cesam and X C les-opal96 - ^cesam respectively. 



Case M/Mq Type AR/R AL/ L Ap c / p c AT C /T C 









A 


B 


C 


A 


B 


C 


A 


B 


C 


A 


B 


C 


Cl.l 


0.9 


MS 


-0.01 


0.10 


-0.08 


-0.03 


-0.05 


-0.02 


0.00 


-0.17 


-0.17 


-0.01 


-0.02 


-0.007 


C1.2 


1.2 


ZAMS 


-0.05 


-0.12 


-0.07 


0.04 


-0.06 


-0.10 


-0.001 


-0.17 


-0.17 


0.006 


-0.04 


-0.05 


C1.3 


1.2 


SGB 


0.01 


0.33 


0.31 


0.36 


0.90 


0.53 


-0.30 


-2.5 


-2.2 


-0.02 


0.55 


0.57 


C1.4 


2.0 


PMS 


-0.16 


-0.13 


0.04 


0.26 


0.20 


0.04 


-0.16 


0.47 


-0.30 


0.00 


0.00 


0.00 


C1.5 


2.0 


TAMS 


-0.03 


0.43 


0.46 


0.76 


0.82 


0.07 


-0.20 


-0.30 


-0.10 


0.06 


0.05 


-0.007 


C1.6 


3.0 


ZAMS 


-0.14 


-0.14 


0.00 


0.14 


-0.14 


0.00 


-0.06 


-0.05 


0.00 


0.01 


-0.008 


0.01 


C1.7 


5.0 


MS 


0.05 


0.18 


0.13 


0.003 


0.25 


0.25 


-0.005 


-0.05 


-0.05 


0.000 


0.02 


0.02 
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Fig. 8 Plots in terms of the relative radius of the differences at fixed relative mass (two left panels) for the internal regions, and at fixed radius 
(two right panels) for the outer layers, between models computed with different opacity routines for the Cases 1.1, 1.2 and 1.3. Solid lines 
correspond to the difference between two types of CLES models: those obtained with the standard CLES version that uses the smoother OPAL 
opacity tables (opal96-S) and those obtained by using an OPAL opacity table obtained without smoothing. Dotted lines: differences between 
the standard CLES models and those from CESAM. Dashed lines: differences between CLES models computed by using OPAL opacity tables 
without smoothing and CESAM models. Left panel: logarithmic sound speed differences. Central left panel: logarithmic pressure differences. 
Central right panel: logarithmic sound speed differences. Right panel: logarithmic adiabatic exponent differences. 



ing to the same radius an oscillatory component (amplitude 
0.01 liHz) remains in the Av because of the differences in 
i~i in the outer layers. 



4 Nuclear reaction rates 

We used the basic pp and CNO reaction networks up to 
the 17 0(p,a) 14 N reaction. In the present models the CESAM 
code takes 7 Li, 7 Be and 2 H at equilibrium while CLES fol- 
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Fig. 9 As Fig.[lOl for Cases 1.4, 1.5, 1.6 and 1.7. 



lows entirely the combustion of 7 Li and 2 H. The nuclear re- 
action rates are computed using th e analytical formula e pro- 
vided by the NACRE compilation dAngulo et all 1 19991) . CE- 
SAM uses a pre-computed table0 while CLES uses directly the 
analytical expressions. Comparing the nuclear reaction rates 
for a 2 Mq stellar structure and a given chemical composi- 
tion, we found that the relative differences are of the order 
of 3 x 10 _4 -2 x 10~ 3 (except for nuclear reactions involv- 
ing 7 Li) if screening factors are included, and of the order 
or 1 0~ 8 if no t . In bo th codes weak screening is assumed un- 
der ISalpeteri (Il954lr s formulation. The screening factor is 

written / = exp(AziZ2\J~^) where z\ and Z2 are the charges 



2 By using pre-computed tables the numerical coherence of deriva- 
tives required in the CESAM numerical scheme are guaranteed. 



of the interactin g nuclei. CESAM uses the expression 4-221 
of lClavtonl dl968h where A = 1.88 10 8 , t, = +z,)x„ 

i 

and Xj is the abundance per mole of element i. The stan- 
dard version of CLES code takes A = 1.879 10 8 and E, = 

4 

L ( 1 +Zj)xj +Z(1 +Z)x(Z) where x(Z) is the abundance 

i'=l 

of an "average" element containing all the elements differ- 
ent from hydrogen and helium, and Z is the average charge 
of this element. This approximate estimation of E, has been 
changed in CLES by assuming full ionization and taking the 
contribution from each mixture element into account. With 
this new prescription the differences between CESAM and 
CLES nuclear reaction rates are still of the same order, but 
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Fig. 10 p-mode frequency differences between TASK 1 models produced by CLES with the two kinds of opacity tables (dash-dotted lines: 
smoothed OPAL; and dashed lines: without smoothing) and CESAM. The latter is taken as reference, and the frequencies have been scaled to 
remove the effect of different stellar radii. For each model we plot two curves corresponding to modes with degrees i = and £ = 1. 



with CESAM values larger than CLES ones at variance with 
what was found with the standard CLES formulation. 

All the Task 1 and Task 3 CLES-models were computed 
with the updated version. 



corresponding values of T e g, logg and chemical composi- 
tion. 

While in CLES atmosphere computations the condition 
at the optically thin limit (p(T m ; n )) is determined for each 



5 Atmosphere 

Eddington's grey T(t) law is used for the atmosphere cal- 
culation: T = r e ff[|(T+ where T is the optical depth. 
CESAM integrates the hydrostatic equation in the atmosphere 
starting at the optical depth T = T m i n (T m ; n = 10 -4 for solar 
like models) and makes the connection with the envelope at 
T = 10 where the continuity of the variables and of their first 
derivatives are assured. The radius of the star is taken to be 
the bolometric radius, i.e. the radius at the level where the 
local temperature equals the effective temperature (t = 2/3 
for the Eddington's law). 

In the stellar structure integration CLES gets the external 
boundary conditions ( the values of density and temperature 
at a given optical depth t) by interpolating a pre-computed 
table, and the stellar radius is denned as the level where 
T = T e ff. The Eddington atmosphere table, which provides p 
and T at T = 2/3 (therefore at R = R t ), was built by integrat- 
ing the hydrostatic equilibrium equation in the atmosphere 
starting at an optical depth that can vary between 10~ 4 and 
10~ 2 . The atmosphere structure for a given model is com- 
puted afterwards, by integrating the same equations for the 
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Fig. 11 Stellar structure of the external layers of 2 M & models at two 
different evolutionary stages (X c = 0.50 and 0.01). Dashed lines: CLES 
models. Dotted lines: CESAM models with boundary conditions given 
by P{tmm)- Solid lines: CESAM models with boundary conditions 
given by p(T m j n ) and default values (see the text). 



Grids:CLES+LOSC 




Fig. 12 Plots in terms of the relative radius of the differences at fixed relative mass (three left panels) between three different CLES computations 
and CESAM models. Solid lines correspond to differences between standard CLES and CESAM. Dashed lines correspond to differences between 
CLES models computed by doubling the number of time steps and CESAM ones. Dot-dashed lines correspond to differences between CLES 
models computed with a double number of mesh points. In lower panels there the differences between CESAM and CLESX2 models computed 
using the opacity tables without smoothing (opal96) are plotted by using dotted lines. The right column plots show the distribution of mesh 
points inside the CESAM models (dotted lines) and CLES ones (solid lines). 



(T e ff, logg, X, Z) by a Newton-Raphson iteration algorithm, 
CESAM allows to integrate the atmosphere by fixing either 
p(^min) or P(tmm)- We think it is worth warning here about 
the relevance of an appropriate choice of p(T m i n ) in CESAM 
calculations. As explicitly indicated in the corresponding tu- 
torial, the default values were determined for solar like mod- 
els, and if much different physical conditions are considered, 
the boundary conditions in the optically thin limit should be 
coherently changed. 

These comparisons have allowed to show the discrepan- 
cies in frequencies that a "black-box" use of an evolution 
code might lead to. For illustration, in Fig.[TT]we show the 
outer layers of 2 M© models at two different evolutionary 
states (X c = 0.50, and X c = 0.01), with effective temperature 
r ef f=8337 K and 6706 K respectively. While the atmosphere 
and sub-phostospheric structure of CLES and CESAM models 
is quite close if P(T m i n ) option is used in CESAM, discrepan- 
cies that increase with the effective temperature of the model 



appear when p(T m ; n ) with default values (derived for solar- 
like models) is adopted. The differences induced in the struc- 
ture of these outer layers by the use of inappropriate limit 
values in the atmosphere integration are much larger than 
those due to opacity differences at logT ~ 4. These outer 
structure differences can lead in fact to frequency differences 
of the order of several fiHz. 



6 Numerical aspects 

The different numerical techniques in CLES and CESAM lead 
to different distribution of mesh points in the stellar struc- 
ture and to different values of the time step between two 
consecutive models. A s pointed out in Sect. 13.11 (see also 
iLebreton et al. ., 2008a) the disagreement between CLES and 
CESAM models can be partially reduced in some cases by 
changing the mesh. Even if both sets of models have a sim- 
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ilar total number of mesh points, their distribution inside 
the star, as shown in Fig. [T2] (right column), is quite dif- 
ferent. In this section we analyze the effect of doubling the 
number of mesh points (CLES-X2) or the number of time 
steps (CLES-T2) on the differences between CLES and CE- 
SAM. As shown in Fig. [12] the effects of these changes are 
not the same in all the considered stellar cases. While for 
the CASE 1 . 1 the increase of mesh points makes almost dis- 
appear the disagreement between CLES and CESAM models 
(51nVar ~ 10~ 4 — 5 x 10~ 5 ), the effect is almost negligible 
for the CASE 1.3. Doubling the number of mesh points in 
CASE 1 .5 leads to a significant effect in decreasing the lu- 
minosity of this TAMS model. As pointed out in Sect. 13.11 
a larger number of mesh points near the boundary of the 
convective core decreases the effect of the sort of "numeri- 
cal diffusion" that changes the chemical composition gradi- 
ent at the boundary of the convective core, and that works 
as a slightly larger overshooting. In fact, the differences of 
hydrogen abundance in the region of chemical composition 
gradient (r/R between 0.035 and 0.06, Fig. 1121 central-left 
panel for CASE 1.5) also decrease with respect to those ob- 
tained with the standard CLES models. As already discussed 
in Sect. 13.11 a part of the disagreement between CLES and 
CESAM comes from the differences in the opacities used in 
both codes. In the lower panels of Fig.Q~2]we have also plot- 
ted (dotted lines) the results of comparing the models com- 
puted with CLES doubling the number of mesh points and 
using the OPAL96 tables without smoothing. A significant 
decrease of luminosity and hydrogen-profiles differences is 
obtained when both the number of mesh points and the opac- 
ity tables are changed. 

Decreasing the time step in the evolution models does 
not lead, in general, to better agreement between CLES and 
CESAM. 



7 Conclusions 

In addition to the quantitative re sults of code comparison 
presented in lMonteiro et al\ d2008l) and lLebreton et al.\ d2008ah 
the analysis of stellar models computed with the codes CE- 
SAM and CLES has allowed us to reveal some interesting as- 
pects about the two codes, as well as about the input physics, 
that only a thorough analysis might bring to light. Some of 
these evidences have led to changes or correction of bugs in 
the codes, other simply allowed us to understand the origin 
of some differences that were reported in the above men- 
tioned papers. 

- The inconsistencies among the thermodynamic quanti- 
ties in OPAL2001 equation of state tables lead to dif- 
ferences in the stellar models and in the oscillation fre- 
quencies larger than the uncertainties due to different in- 
terpolation tools. Even if the quantity Cy as tabulated 
in OPAL2001 tables is not used, the remaining incon- 
sistencies among the three adiabatic indices lead to dif- 
ferences between the model computed with a code that 
takes the thermodynamic quantities directly from OPAL 



tables (such as CESAM) and a model computed with a 
code whose thermodynamic variables are derived by the 
thermodynamic relations and a minimum of tabulated 
quantities (such as CLES). Furthermore, the discrepan- 
cies will depend on the choice of tabulated quantities. 

- The precision of the theoretical oscillation frequencies 
is seriously limited by the uncertainties in the opacity 
computations. 

- Different approaches used to estimate the electron den- 
sity in CESAM and CLES lead to differences in the screen- 
ing factors that have no relevant effects on the stellar 
models. 

- Even with a simple physics, such as the Eddington's law 
for gray atmosphere, the details of numerical tools can 
have significant consequences on the seismic properties 
of the models. 

- The different distribution of mesh points in the models 
can explain part of the disagreement between CESAM 
and CLES models. An increase of mesh points in the in- 
ternal regions seems to be required in CLES to decrease 
the differences with CESAM. 

Apart from the discrepancies in the screening factors which 
does not significantly affect the oscillation frequencies, the 
other factors analyzed here can affect the absolute oscillation 
frequencies by up to several /iHz. 
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